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Abstract 

The mass spectra of heavy and light mesons is computed within the framework of the relativistic 
flux tube model. A good agreement with the experimental data is obtained provided that the flux 
tube contributions, including retardation and spin-orbit effects, are supplemented by a one-gluon- 
exchange potential, a quark self-energy term and instanton-induced interactions. No arbitrary 
constant is needed to fit the absolute scale of the mass spectra, and the different parameters are 
fitted on lattice QCD in order to strongly restrict the arbitrariness of our model. The relevance of 
the present approach is discussed in the case of glueballs, and the glueball spectrum we compute 
is compared to the lattice QCD one. Finally, we make connections between the results of our 
model and the nature of some newly discovered experimental states such as the /o(1810), X(3940), 
Y(3940), etc. 
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I. INTRODUCTION 



Potential quark models have been proved to successfully reproduce the experimental 



Si 



ince the 



meson and baryon mass spectra for nearly thirty years [1 
pioneering works on this subject , they have been a matter of constant interest. In 
particular, new hadrons are continuously being discovered, and potential models offer an 
intuitive and efficient way to understand the physical properties of these new experimental 
states. Apart from many other relevant effective approaches of QCD Q], lattice QCD has 
recently emerged as a powerful method to deal with the full QCD theory (see Ref. 9] for 
an introduction). Interestingly, some basic features of potential models in the heavy meson 
sector have been confirmed by lattice QCD calculations. The potential energy between a 



static quark and antiquark has indeed been shown to be roughly compatible with a 
of the Cornell form ar — n/r + C id Illfl. which is widely used in potential models 



potential 



12|. The 



total interaction is then seen as the sum of a long-range part which encodes the confinement 
and a Coulomb term arising from one-gluon-exchange process between the quark and the 
antiquark. Finally, C is an arbitrary negative constant used to fit the absolute scale of a 
given mass spectrum. 

The structure of light meson spectra can also be understood by using potential models {2} • 
In particular, a semirelativistic quark kinetic term of the form \/p 2 + m 2 allows to deal 
with li ght , and even massless, particles. As suggested by the background perturbation 
theory [131 ]. the relativistic, spin- dependent, corrections, can then be developed in powers of 
l//i 2 , ft being defined by (^Jp 2 + m 2 ^. In this expression, the average values are typically 
computed with the eigenstates of a spinless Salpeter Hamiltonian with the Cornell potential. 

u n 

Such a framework leads to interesting results concerning light mesons [7|, H4|], and reduces 
to the usual formalism for heavy mesons, since fi ~ m at large quark mass. Developing 
the relativistic corrections in powers of l//i 2 is rather natural: It is actually nothing but 
developing them in terms of the quark energies rather than in terms of their masses. Such a 
procedure has already been successfully applied to light mesons in many previous potential 
models of mesons and baryons 

A more recent approach, called the relativistic flux tube model, is an effective meson 



model obtained from the full QCD theory 



15 



181 ] . It relies on the assumption that 



the quark and the antiquark in a meson are linked by a straight color flux tube, actually 



2 



a Nambu-Goto string, carrying both energy and angular momentum. This string, or flux 
tube, is responsible for the confinement. It is worth noting that lattice QCD simulations 
show that the chromoelectric field between a static quark-antiquark pair is roughly constant 
on a straight line joining these two particles 191 ]. validating the physical picture of the 
relativistic flux tube model. Like potential models, it produces linear Regge trajectories in 
the ultrarelativistic limit, that is a linear link between the square masses of light mesons 
and their angular momentum, and it reduces to the usual Schrodinger equation with a linear 
potential in the heavy quark limit. We made in previous works several attempts to include 
in the relativistic flux tube model physical mechanisms that were usually neglected: Quark 



self-energy 



retardation effects 



21 



22 ]. and spin interactions [231 ] . Apart from these 



corrections, that are all due to the flux tube itself, it is also necessary to include short-range 
potentials such as the Fermi-Breit interaction, which comes from the one-gluon-exchange 
diagram between a quark and an antiquark, in order to get a satisfactory description of 



mesons with the relativistic flux tube model 



24 



25]. 



Our aim in the present paper is to put together all these physical contributions, which 
were up to now only separately studied. By doing this, we will build a generalized relativistic 
flux tube model leading to mass spectra in quantitative agreement with the experimental 
data, and to relevant suggestions concerning the nature of some recently discovered states 
such as the / (1810), X(3940), F(3940), ...This model is presented in Sec. M Then, 
we apply it to compute heavy and light meson mass spectra in Sees. IIHI and IIVI where 
we systematically compare our data to the available experimental states. The application 
of potential models t o g lueballs - bound states of gluons - has already been done in the 



past [26|, [27|, [28|, |29j, [30), |3l|, in good agreement with lattice QCD results. However, the 
validity of such models to describe glueballs in still controversial. We discuss in Sec. |V] the 
possible extension of our relativistic flux tube model to the case of a bound state of two 
gluons. Finally, we study the mixing between scalar mesons and glueball in Sec. I VI I and 
draw some conclusions in Sec. IVHI 
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II. THE MODEL 



A. Hamiltonian 

Starting from the full QCD theory, a Lagrangian for a system of two confined spinless color 
sources can be derived: The corresponding model has been called the rotating string model 
[15I ]. In this approach, originally introduced as an effective meson model, the quark and 
the antiquark are linked by a Nambu-Goto string, or flux tube. This object, characterized 
by its energy density a, simulates the exchange of gluons responsible for the long-range 
part of the interaction, that is the confinement. The rotating string model is completely 
equivalent to another phenomenological description of mesons, the relativistic flux tube 
model . once the auxiliary fields (or einbein fields) appearing in the rotating string model 



are properly eliminated 



13, 



181 ] . Although it can be numerically achieved, the resolution of 



the nonlinear coupled equations of the relativistic flux tube model is a difficult problem [32 ]. 
More convenient expressions can be obtained by using an accurate approximation of this 
approach, that we previously called the perturbative flux tube model The present 

section is devoted to the presentation of the particular Hamiltonian underlying this model. 

The complexity of the relativistic flux tube equations comes from the fact that the contri- 
bution of the flux tube does not reduce to a simple static potential: This object is dynamical, 
and carries orbital angular momentum as well as energy. When the orbital angular momen- 
tum of the meson, denoted as £, is equal to zero, the two-body relativistic flux tube model 
becomes a spinless Salpeter Hamiltonian with a linear confining potential ar. Such a Hamil- 
tonian has been intensively used in the literature to study the properties of hadrons (see a 
review in Ref. [jjj]). As a is the energy density of the flux tube and r is the quark- antiquark 
separation, it is readily observed that the linear potential is actually the total energy of a 
static straight string whose length is r. The dynamical contribution of the flux tube is actu- 
ally small enough to be treated as a perturbation of the spinless Salpeter Hamiltonian with 

nnn 

potential ar, and one then obtains the so-called perturbative flux tube model [14J, [201, 1301 ] . 
If £ is not too large (typically £ < 6), this a ppr oach reproduces the "exact" flux tube mass 



spectrum with an accuracy better than 5% 



20] 



For a bound state of two particles of the same current mass m, that is the case we restrict 
to in this paper, the perturbative flux tube model is defined by the following spinless Salpeter 
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Hamiltonian 

H = 2 v /p 2 + m 2 + ar, (1) 

completed by a perturbation which encodes the dynamical contribution of the string through 
its angular momentum [^(J 

AHur = _ <<+D (2) 
fxr(6jj, + ar) 

The quantity fi appearing in the above equation can be seen as a dynamical quark mass 
given by 

fi=(y/p 2 + m*), (3) 

in which the average value is computed with an eigenstate of the Hamiltonian ([T]). The 
dynamical mass is then state dependent. The contribution AM^ to the total mass can be 
computed with a good accuracy by the following approximation Q 

atU+l) /1\ 
fj,{6fji + a{r)) \r / 

Equations (HI) and O) define the original perturbative flux tube model, as studied in 

n □ 

Refs. [14J, 120( . However, important improvements can be brought by dropping the various 
simplifying assumptions underlying this approach. Three hypothesis are actually used to 
establish the equations of the rotating string model, and equivalently of the relativistic flux 



tube model 18j. The first one is the assumption that the flux tube is a straight line linking 
the quark and the antiquark. As we already pointed out, such an ansatz is in agreement 
with lattice QCD simulations [19]. Moreover, other calculations within the framework of 
effective models show that the possible deviations of the string from the straight line are 
negligible 3, 33] in mesons. 

It is worth mentioning that we only focus on usual mesons in the present work. The flux 
tube model can actually succesfully describe hybrid mesons as well, i.e. mesons in which 
the gluonic field is in an excited state. But, in the case of hybrid mesons, the linear flux 
tube picture does not hold anymore: An excitation of the gluonic field deeply changes the 
configuration of the system. A first possibility of modelizing it within the flux tube model 
is to treat the fluctuations of the flux tube in a quantized way. Then, the different energy 
levels of the flux tube define configurations corresponding to hybrid mesons j^]. A second 
way is to assume that the excitation of the gluonic field creates a constituent gluon [35 ]. 
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A hybrid meson is then seen as a quark- ant iquark-gluon bound state in which the gluon 
is linked to the quark and to the antiquark by two straight strings. In every case, we can 
conclude that a - weakly deformed or not - straight line configuration of the flux tube is 
only valid for usual mesons, when the gluonic field is in its ground state. This is the case 
that we will treat in the latter. 

The second simplification which is made to obtain the relativistic flux tube model is 
the neglect of the retardation effects. It can be shown that these effects can be computed 



in perturbation within the framework of the rotating string model [2l|, by allowing the 
quark and the antiquark to have a different temporal coordinate. Retardation effects bring 
a contribution to the total mass which reads 

AM ret = (5a) 
la?. 



with 



' \12 + fir 2r 



a 2 ap 2 . a \i 



90 20 fir 12 r 

ary 

2 12 + ~40 



y 2 , (5b) 



r{/j + ar/b) 

In these last formulas, p r is the quark radial momentum, while y is its transverse speed. The 
various symmetrizations on the noncommuting operators have not been written in order to 



21] for a detailed study of the retardation 



clarify the notations. We refer the reader to Ref. 
effects, including the computation of the relations (j5J). However, for our purpose, it is 
interesting to notice that retardation brings a negative contribution to the total mass. The 



physical content of Eq. (joajl appears more clearly in two limits: 

3a . , I a \ 



AM ret \ m=Q ^--, AM^I^oc^-J . (6) 

As expected, the retardation effects vanish for heavy quarks, where the typical speed scale 
is very small with respect to the speed of light. On the contrary, they are maximal for light 



particles. An interesting feature of AM ret \ m=Q is that it preserves the Regge trajectories [21 ]. 

The third simplification is that the rotating string model neglects the spin of the quark and 
the antiquark. An attempt to exactly include the spin degrees of freedom in the relativistic 



flux tube model has been made in Ref. 



361 ] . but it leads to very complicated equations, 



that, to our knowledge, have not been solved yet. However, the spin contribution can be 
described in perturbation by a spin-orbit coupling between the angular momentum of the 
flux tube and the spin of the confined particles. In a two-body system where the constituent 
particles have the same mass, this term reads [23 ] 



AH a 



aL ■ S 

2jj?r 



AM, 



L-S 



(?) 



where S = S± + S2 is the total spin of the system, and where L is the relative orbital angular 
momentum in center of mass frame. This formula agrees with the results of back grou nd 
perturbation theory and with those obtained with the Wilson loop technique 10] in 
the limit of heavy quarks. Such a spin-orbit coupling can actually be thought as a Thomas 
precession of the spin of the confined particles in the color field. 

Apart from these flux tube corrections, it has recently been shown that the quark self- 
energy contribution, which is created by the color magnetic moment of the quark propagating 
through the vacuum background field, adds a negative constant to the hadron masses 37J . In 
the case of a meson made of two quarks of the same mass, the quark self-energy contribution 
to the meson mass is given by [20, Q 

/ ar)(m/S) 



AM, 



que 



7T/i 



(8) 



The 77-function is such that 77(0) = 1, and its value decreases monotonically toward with 
increasing quark mass (its explicit form can be found in Ref. 37|). 5 is the inverse of the 
gluonic correlation length, and its value is estimated at about 1.0 — 1.3 GeV by lattice 
QCD computations 38|, [39]. As the quark self-energy contribution varies very little with 
this parameter (less than 1%), its value can be fixed at 1 GeV 2(3]. The factor f has been 



computed by lattice calculations. First quenched calculations gave / = 4 38(. A more 



recent unquenched work 



391 ] gives / = 3, the value that we choose in this work. 



Finally, a significant contribution to hadron masses is also given by the one-gluon- 
exchange mechanism between color sources. These are short-range interactions, which are 
consequently not encoded in the relativistic flux tube model. The one-gluon-exchange po- 
tential for mesons is equal to the well-known Fermi-Breit interaction, that is {40], p. 239], 
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vr^i « 2 -3^ >, ( 9) 

where the Darwin term has been neglected, as is a phenomenological strong coupling 
constant, and U(r) is the gluon propagator in position space. For massless exchanged 
gluons, one simply has U(r) = 1/r, but Eq. (Q remains valid for any form of the gluon 



propagator 



4l|. 



By putting together all these ingredients, a mass spectrum can be numerically computed. 
The numerical method that we use in this work is the Lagrange-mesh method, a remarkably 
simple and accurate method to solve nonrelativistic as well as semirelativistic eigenequa- 
tions 42J. The procedure is the following. Firstly, we compute the eigenvalues Mo of the 
Hamiltonian 



H = 2^/p 2 + m 2 + ar - -a s U{r), (10) 

o 

which is basically a spinless Salpeter Hamiltonian with a Cornell potential, encoding the 
dominant contributions of our model. The ar term comes from confinement, while the 
— (4/3)acsU(r) potential is the lowest order contribution of the one-gluon-exchange poten- 
tial (EJ). Secondly, the dynamical mass /i can be computed thanks to the relation ([31). 
Thirdly, the various contributions we presented here can be added in perturbation to get 
the total mass M, defined by 

M = M + AM str + AM ret + AM so + AM qse + (v oge {r) + ^a s U(r)^ . (11) 

Clearly, M depends on the radial and orbital quantum numbers n and £, but also on the 
spin degrees of freedom. The spectroscopic notation (n + l) 2S+l £j will be used in order 
to unambiguously refer to a particular state. We recall that the parity and the charge 
conjugation of a meson are given by P = (— and C = (— l) e+s . 

Let us stress two interesting advantages of the flux tube model we presented in this section. 
Firstly, it clearly appears that the mass formula we obtain is valid for any value of m. The 
reason is that the different corrections to the dominant Cornell potential are developed in 
powers of rather than 1/m. For heavy quarks, fi ~ m, and we recover usual relativistic 
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corrections of order 1/m 2 . However, there is a hope that this model can be generalized for 
light quarks since, even if m = 0, \i remains approximately equal to 0.3 GeV [30]. It is worth 
mentioning that older nonrelativistic quark models phenomenologically used a dynamical 
mass around this value of /i, with a satisfactory agreement with experiment [11] . Secondly, 
we do not need to add a negative constant to the Cornell potential in order to fit the abso. 
mass scale of the spectrum, as it is the case in most of the potential models 
Formally, this constant can be thought to be replaced by the recently found quark self-energy 
and retardation terms, which are state-dependent. 



ute 



12|. 



B. Effective quark size 



Within the framework of a potential QCD model, it is natural to assume that a quark is 
not a pure point-like particle, but an effective degree of freedom that is dressed by a gluon 
and quark-antiquark pair cloud. Such an hypothesis leads to very good results in the meson 
4\ and baryon ^] sectors, and can also be applied to glueballs 29, 43]. As in these two last 
references, we assume a Yukawa color charge density for the confined particles 



p{u) 



-u/y 



47T7 2 



u 



(12) 



where 7 is the effective size. It is readily checked that p(u) reduce to 5 3 (u) when 7 tends to 
zero. The interactions between the two confined particles are then modified by this density, 
a bare potential V(r) being transformed into a dressed potential V(r). This potential is 
obtained by a double convolution over the densities of each interacting particle and the bare 



potential. It can be shown that, when both confined partic 
this procedure is equivalent to the following calculation 



29 



es have an equal effective size 7, 



44 1 



V(r) 



d 3 r'V{f')T(f-f'), withT(u) 



-u/y 



57T7 



(13) 



This is actually our case, since we only deal with bound states of two particles of the same 
mass. So, we can assume that their effective sizes are identical. Other color-charge densities 
could be used, a Gauss la n one for instance fl. We have nevertheless st r o„ g indications 
that such a change cannot noticeably modify the results 45[ ; moreover, all convolutions are 
analytical with this Yukawa form. 

The application of the convolution to the Fermi-Breit potential ([9]) can be achieved by 
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replacing the gluon propagator U(r) by a "dressed" one, that is the convoluted propaga- 



tor 



U{r) 



-r/l 



(14) 



1 1 

r 2-i / 

This completely removes all singularities in the short-range interactions. For consistency, 
we apply the same transformation to the confining terms, through the substitution [29| 

4a7 2 



ar 



ar = ar + 



-r/7\ _ 



) — a^e 



(15) 



C. The strong coupling constant 

Lattice QCD clearly shows that the static potential between a quark and an antiquar 



k is 

compatible with the Cornell form ar — (4/3)asr, for a ~ 0.20 GeV 2 and as ~ 0.22 l_l|, p. 
42]. Consequently, the total energy is separated into a confining part, and a "residual" short- 
range part - let us note that this separation is de facto performed in our model. The Coulomb 
term is clearly equal to the lowest order approximation of the Fermi-Breit potential, with 
a rather small value for as- Thus, the effective strong coupling constant does not blow up, 
even in the bound states sector we are dealing with, provided that nonperturbative effects 
are correctly taken into account. Several attempts have been made in order to include 
the nonperturbative contributions into the well-known formula giving as{q 2 ), q being the 



transfered 4-momentum 



46j. All these approaches qualitatively lead to: 



2 127T 

(33 - 2N f ) In (^±p^) 

where £(<? 2 ) is a monotonic function such that £ 2 (0) > and £ 2 (<7 2 — > oo) — > 0. Nf is the 
number of quark flavors whose masses are lower than q 2 , and A is the famous lambda QCD 
parameter. However, the explicit formula giving £(q 2 ) is different following the different 
works. Equation (fl6l) states the strong coupling constant remains finite for q 2 = A 2 and 
tends to a maximal value «s(0) < 1, in agreement with lattice QCD j^ . 

A simple, phenomenological, way to mimic the behavior of Eq. ffTUj) in position space is 



to replace as by [47] 



a 5 (r) = « (l-e-^). (17) 

It is nowadays well-established from experimental measurements that as{q 2 = rn 2 z ) = 
0.1176 ± 0.0020 [48], with rriz ~ 91.19 GeV. The parameter r c can consequently be fixed by 
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demanding that as{q 2 



rrir 



a sir = \jvnz). This condition leads to 



mz In 



a s {q = rrq 
a 



-i 



(18) 



The only remaining parameter is cto, that we will fix by a comparison of our results to lattice 
QCD in the following section. 



D. Summary of the model 

We have now introduced all the ingredients of our model, and one may wonder about 
the possible double counting of the interactions. Let us sum up the different contributions 
which are taken into account. 

First of all, the long-range (confining) interactions are encoded in the flux tube. At the 
lowest order, it reduces to a linear potential, and is thus independent of the features of the 
confined quarks. As the flux tube is a dynamical object, it carries angular momentum; its 
dynamical contribution to the total mass, AM r f t , can be added perturbatively, so defining 
the perturbative flux tube model. Moreover, the flux tube is assumed to emerge from a flux 
of exchanged gluons, and retardation effects due to the finite value of light speed are thus 
present throu gh t he term AM ret . Then, the spin of the quarks comes into play. It has been 
shown in Ref. 37J that the interactions between the spin of a quark and the flux tube can be 
split in two parts. The first one involves the dynamics of the flux tube through its angular 
momentum, and leads to the spin-orbit term AM so . This term is a relativistic correction 
that can be thought as a Thomas precession in the color field. The second part does not 
depend on the flux tube momentum and, after tedious calculations, give rise to the quark 
self-energy contribution AM qse [371 . Quark self-energy and spin-orbit are thus two different 
aspects of the interactions induced between the spin of the confined quarks and the flux 
tube, and both contributions are distinct from dynamical and retardation terms, which are 
spin-independent. 

As the flux tube model completely neglects the short-range interactions, these have to be 
added as a Fermi-Breit potential, V oge (r), coming from one gluon exchange diagram. In the 
same way, the self-energy diagram, which, among other QCD Feynman diagrams, causes the 
running of the coupling constant, is not contained in the previously mentioned quark self- 
energy term, which is only linked to the flux tube. Short-range interactions thus complete 
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the model by adding interactions which are all neglected by the flux tube. Finally, mass 
formula (fTTj) will be used to compute all the mass spectra in the following. But, V(r), as, 
and U(r) will be replaced by V oge (r), as(r) and U(r) in this last formula, in order to take 
into account the effective quark size and the running of the coupling constant. 



III. HEAVY MESONS 



A. Fitting the parameters 



Potential quark models for heavy mesons have proved to be particularly successful for 
a long time, and still deserve interest because of the new heavy states which are currently 



being discovered [2|, |6|, |49j, |50| . Moreover, even the complete Fermi-Breit potential, including 
the relativistic corrections, has recently been validated by lattice QCD calculations 5l|. The 
various levels of the potential energy between two static quarks have also been computed in 
521 ] . As we already pointed out, the ground state level, denoted as the E+ one 



lattice QCD 



in Ref. 



521 ]. appears to be compatible with a standard Cornell potential. But, in order to 
consistently fit our parameters, we actually have to fit this ground state with the following 
convoluted potential 



V h (r) =ar - -a s (r)U(r). 



(19) 



It is actually the potential coming from our complete model, where the relativistic corrections 
vanish in the static quark limit (m — > oo) since they are expressed in powers of l//x ~ 1/m. 
As we can see in Fig. [T], the best fit is obtained with the standard values a = 0.185 GeV 2 , 
a = 0.400, and jh = 0.200 GeV -1 (the h index denotes the heavy quarks). Such a value for 
«o is rather usual in potential models, and defines the value r c = 0.031 GeV^ 1 through the 
relation (fTBl) . The parameters of potential (TT9T) are thus completely fixed by lattice QCD. 
Let us note that a is usually assumed to be around 0.19 GeV 2 [ll, p. 9], in agreement with 
our fit. 

For heavy quarks beyond the static limit, relativistic corrections come into play. They do 
not demand any new parameter since they are expressed in powers of l//z, given by Eq. ([3]). 
However, two parameters must be added because of the quark self-energy. Following the 
latest lattice QCD results, we set / = 3.000 and S = 1.000 GeV j^l, although this last 
parameter has a very little influence (less than 1%) 0]. Finally, the masses of the c and 
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b quarks will be fitted in order to reproduce the J/ip and T(1S) mesons respectively. The 
values of the different parameters which are used are summarized in Table HI 

Before computing detailed mass spectra, it is interesting to make some comments about 
the influence of the different terms appearing in the mass formula (llip . The numerical 
values of the different contributions in this last relation are given in Table [III for several 
qq states, using the parameters of Table [H As it can be observed, the retardation term, 
AM ret , always brings a negative contribution around 0.1 — 0.2 GeV, which is maximal for 
light mesons and for small quantum numbers. The parameters which have to be used in 
the case of light mesons will be discussed in Sec. IIV} but performing some calculations 
are already interesting at this stage in order to understand the evolution of the relativistic 
corrections when one deals with light instead of heavy quarks. The dynamical contribution 
of the relativistic flux tube, AM r f t , also decreases the mass of the mesons, and is logically 
more and more relevant at high £, since the flux tube carries orbital angular momentum. A 
last term decreases the meson mass in every case: The quark self-energy one, denoted as 
AM qse . Its contribution is very important for the nn states (up to 0.5 GeV), but quickly 
decreases for heavier mesons. We recall that the symbol n is used to denote both light u or d 
quarks. It becomes negligible for bb states. The spin-orbit term coming from the flux tube, 
AM so , does not contribute significantly for heavy mesons, but can become relevant for nn 
states. Finally, the relativistic, short-range, corrections ^given by (v oge (r) + ^as(r)U(r)X) 
bring state-dependent contributions, ensuring the mass splitting of the various J PC states. 



B. Mass spectra 

We begin our study by computing the bb meson spectrum, corresponding to the T meson 
family. As it can be observed in Fig. [2, the agreement between the experimental data and 
our spectrum is rather satisfactory, for a b mass given by m?, = 4.785 GeV. The radial 
trajectory of the T mesons is well reproduced by computing the 1 , (n + l) 3 Si, states . 
Moreover, our results are compatible with the interpretation of the Xbj mesons as the l 3 Pj 
and 2 3 Pj triplets, although we find them slightly lower than the experimental data. 

Among the heavy mesons, charmonia have recently become an important source of discov- 



eries of new experimental states such as the V(3940), V(3940), etc. (see Refs. 



53 



54 



55|). 



That is why a computation of the cc spectrum within the framework of our generalized flux 
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tube model is of particular interest. The results are plotted in Fig. [HJ It appears that our 
spectrum fits the data with a reasonable agreement, excepted for the ■0(4040) and ■0(4415) 
mesons, which are expected to be a 3 3 5*i and a 4 3 S*i state respectively. They are overes- 
timated in our model, as it is the case in a previous study relying on a Cornell potential, 
thus sharing some similarities with this work 56]. The ^(3770) is however well matched by 
a 1 3 D% state. We can also observe that the predictions of our model concerning the 2 3 Pj 
cc multiplet, for which no firm experimental candidate is known, are compatible with the 
lattice QCD study of Re f. |57l| . as well as with other studies on the charmonium spectrum 
in lattice QCD (see Ref. [581] for a review on the subject). 

It is now interesting to focus on the recently discovered X(3940), F(3940), and Z(3930) 
states. It has been argued in the literature that the V(3940) and the Z(3930) are the 
Xci(2P) [50] and the Xc2(2-P) 55] respectively. Our results agree with that suggestion. But, 
it has also been sues; ested that the X(3940) is the Vc (3S) 0, fl. However, as it has 
already been pointed out, quark models predict a too high mass for that state 55[. The 
most satisfactory assignment for the X(3940) from the point of view of our flux tube model 
is that this state is the Xco(2-P), in agreement with our spectrum (see Fig. [3]). Consequently, 
we suggest that these three X, Y, Z states could be identified with the x c j(2P) triplet. 



IV. LIGHT MESONS 

The model we presented in Sec. [TT] has the nice feature that it is well-defined for light 
particles, even for massless ones. That is why the computation of light meson spectra within 
our framework seems to be relevant. 



A. Fitting the parameters 

It is worth discussing a bit about how the parameters of the model should be modified 
in the light meson sector. Rigorously, an effective model of light mesons should take into 
account the possible decays and coupling with other channels for any given state. This is 
clearly not the case in our model, as in most of the quark models of mesons. Generally, this 
problem is bypassed by fitting the parameters on the experimental spectrum in order to have 
an agreement as satisfactory as possible between theory and experiment. Up to now, we 
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did not fit the parameters on the experimental data, but rather on lattice QCD. Actually, 
our way of proceeding is indirectly based on experiment through lattice QCD. It is indeed 
known that this last approach correctly reproduces the heavy meson spectra (see for example 
Ref. [581] ). But, its advantage is that it can also provide us with an estimation of the quark- 
antiquark potential energy and of the values of / and 5. It is thus more straightforward to 
fit our potential on the corresponding lattice QCD curve than on many experimental data, 
knowing that lattice QCD correctly reproduces the experimental spectra. By this way, we 
are led to particular values for a, ao, and 7/,, but the quark masses m;, and m c have still to 
be adjusted on the experimental masses, as it is usually done in potential models. 

Neither the flux tube energy density a nor the and the coupling constant at zero momen- 
tum, ao, are expected to vary with the quark mass, thus the same values than in the case 
of heavy mesons can be kept to describe light mesons. The same remark holds for / and S. 

We turn now our attention to the light quark masses. The mass of the n quark, is fixed 
at m n = for future convenience, in very good approximation of the current experimental 
data, stating that 1.5 MeV< m u ,md < 7.0 MeV 48J. Let us note again that no distinction 
has to be made between u and d quarks because our model is isospin-independent. The 
mass of the s quark is fitted to reproduce the 0(1020), which is largely accepted as a pure 
ss state (see for example Ref. [2}). We obtain m s = 0.202 GeV. 

Finally, we have to ask whether the effective size 7^ should be modified or not in the 
light quark sector. Intuitively, one can imagine that a heavy quark is closer of the classical 
picture of a pointlike particle than a light one. That is why it seems justified to assume that 
li > Jh, where the I index stands for light quarks. After a fit on the experimental data, we 
find indeed that the best agreement with our model is achieved for 7; = 0.940 GeV" 1 . That 
is thus the value that we will choose. Interestingly, such a value corresponds to the effective 
quark size which is generated by instanton-induced interactions, as we will mention in the 
next section. 



B. Instanton-induced forces 

The lightest mesons are the nn ones. Since we set m n = 0, it is worth mentioning that 
the results are exactly independent of the value of 5. An additional physical mechanism 
has to be taken into account in the nn mesons, which is the instanton-induced forces. 
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Instantons are classical solutions of the euclidean equations of motion of QCD, which provide 
informations on the nontrivial vacuum structure in QCD. In a light meson, it has been 



shown that instantons induce forces between the quark and the antiquark (see Ref. 59( for a 
review) . Such forces can be included as a potential in quark models, and only give a nonzero 
contribution for the 1 5'o nn states j^j]. A proper calculation of the influence of these effects 
on the light meson mass spectra is rather complicated, and was performed in Ref. j^] with a 
set of parameter nearly identical to ours, denoted as Model III. The main change resides in 
the value of the string tension a. But, if AM ins is the instanton contribution, AM ins / y/a is 
universal since we deal with massless quarks (m n = 0). We can then compute from Ref. 4j 
that 

AM; S = -0.438 GeV, AM^ 300 ' = -0.205 GeV, AAC ( S 1800) = -0.119 GeV, 
AM? ns = -0.042 GeV. (20) 

Instanton-induced forces do not act on other states than the ir and rj ones. They are thus 
very different in nature than the quark self-energy contribution, which acts on every qq 
state, although both terms are related to quark interactions with a nontrivial configuration 
of the gluonic field. Actually, the different status of these interactions could be understood 
by considering the different length scales appearing in the model. Indeed, the flux tube 
configuration emerges at large distances and give rise to numerous interactions, which were 
presented in Sec. Ull including quark self-energy. Instantons are however particular con- 
figurations of the gluonic field whose typical size is well smaller than the typical size of a 



hadron (roughly related to the size of the flux tube) [60(, and they thus lead to different 
phenomenological implications. 

It can be shown that instantons give an effective size to light quarks, which is generally 



39. 



59 



601 ] . This effective size is actually linked to the 



computed to be around 1 GeV -1 
nonzero value of the mixed quark-gluon condensate g (q a^F^q) (see for example Refs. 39 , 
60j). Remarkably, this effective size is close to the fitted value that we use, namely 7; = 0.940 
GeV -1 . The interpretation of our phenomenological effective size in terms of the effective 
size generated by instantons in the light quark sector seems thus rather natural. It is worth 
mentioning that, rigorously, instantons also confer a constituent mass to the light quarks, 
clearly because of the breaking of chiral symmetry. But, this constituent mass (of order 10 
MeV) is small enough with respect to the dynamical mass \i ~ 300 MeV generated by the 

16 
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confinement [4|, and then can be neglected. 
C. Mass spectrum 

The nn spectrum is plotted in Fig. H] and compared with the experimental data. As our 
model neglects the isospin, we only get one point corresponding to each isospin doublet. We 
see that our spectrum is in correct agreement with these data, excepted for the 7r(1800) - 
but the experimental evidence for this last state is weak. It is interesting to notice that the 
mass we find for the lowest ++ state is 1.351 GeV, close to the /o(1370). Our model thus 
support the idea that the / (1370) is a dominantly nn state, as it is argued in Ref. [61 1 . 
Among the numerous 2 ++ mesons, the /2(1950) seems to be a pure nn state, with quantum 
numbers 2 3 P 2 - 

It is worth mentioning that no scalar state is present between 0.616 GeV and 1.346 GeV 
- corresponding to the n meson and to the 7/(1295) and 7r(1300) respectively - in our nn 
spectrum. However, recent experimental results have provided evidence for low-lying states 
with quantum numbers ++ and masses in this range, namely the / (600), k(800), / n (980) 



and a (980) (useful references can be found in [62j). The recent detailed study of Ref. [62 ] 



suggests that these states are not qq mesons, but rather that they are part of the lowest-lying 



scalar tetraquark nonet, in agreement with earlier works on this subject [63|. Following these 
results, we thus conclude that these light experimental scalar states cannot be reproduced 
by our model, since we only deal with qq mesons, and not with qqqq tetraquarks. 

A correct description of ss mesons cannot be achieved without including a mixing with nn 
mesons [2| . As we neglected such couplings, a detailed ss spectrum will not be computed. We 
obtain m s = 0.202 GeV. Then, the lowest ++ state (1 3 P ) has a mass given by 1.528 GeV, 
which is near the /o(1500). This state could then be a mostly ss state, in agreement with 
Ref. [(31]. Moreover, the 2 3 P 2 and 1 3 P 2 states have the following mass respectively: 2.062 
and 2.226 GeV, which suggest that the / 2 (2010) and / 2 (2300) are principally ss mesons. 

Let us end this analysis of light mesons by a computation of some highly excited nn states. 
The interest of performing such a calculation actually comes from future experiments like 
PANDA and GLUEX, which will scan the hadron spectrum in the range 1.5 — 3.0 GeV. 
One can hope that these experiments will not only find new mesons, but also detect more 
intriguing particles like glueballs or hybrid mesons. If we assume that a hybrid meson is 
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a qqg bound state, it can be shown that the lowest-lying ones should have the following 
quantum numbers 64]: J PC = ±+ , 1 ±+ , 1 ±_ , 2 ±+ . Among these quantum numbers, only 
1~ + is an exotic state, that is a state with quantum numbers that a usual meson cannot 
possess. The discovery of such an exotic state could be the unambiguous signature of a 
hybrid meson. On the contrary, the other possible J PC are accessible for usual mesons. It 
seems thus interesting for the interpretation of the future experimental results to compute 
with our model the masses of the ±+ , 1 ++ , 1 ±_ , 2 ±+ nn states in the interval 2 — 3 GeV, in 
order to suggest energy levels around which an usual meson should be located. This could 
eventually lead to identify hybrid mesons with non exotic quantum numbers at PANDA or 
GLUEX. Results are plotted in Fig. [51 Since our model neglects the isospin interactions, 
this plot has to be seen as predictions giving masses around which an experimental nn state 
should be located rather than accurate predictions of meson masses. As an example of a 
possible application of Fig. El we can look at the 1 states. Following our predictions, no 
nn states with J PC = 1 should be present in the interval 2.4 — 2.6 GeV. This defines 
a particular mass range in which it should be interesting to search for a hybrid meson in 
future experiments. Indeed, if a 1 state with isospin 1 was found in this mass range for 
example, it could be a serious candidate for a ring hybrid meson: An isospin equal to 1 
forbids a possible ss state, and no glueball with such J PC is expected at such a low mass, 
as we will see in the following. 



V. GLUEBALLS 
A. Parameters 



The model we present in this study, although being rather simple, can describe with 
a satisfactory agreement the meson spectra, even in the light sector. The next step is to 
wonder whether such a flux tube approach can be applied to glueballs or not. The subject 
is still a matter of controversy. Let us review the different parts of our model, assuming in 
a first time that a glueball can be seen as a bound state of two constituent massless gluons 
(m g = 0). The situation is then very similar to light quarks. The dynamical flux tube and 
retardation terms are only defined by the flux tube itself and do not depend on the nature 
of the confined particles. Moreover, we have recently shown that the spin-orbit term arising 
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from the flux tube has the same form for particles of arbitrary spin 23j. On the contrary, 
the one-gluon-exchange potential should be modified, since the effective potential emerging 
from the Feynman diagrams involving a pair of gluons at tree level is different of the one 
between a quark and an antiquark. It has been computed in Ref. 25[], and reads 



Voge(r[ 
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(21) 



where U(r) is again the gluon propagator in position space, and where n is the dynamical 
gluon mass, which typically has a value around 0.5 — 0.6 GeV, in agreement with other 



effective approaches 30, [65|]. A last modification occurs in the self-energy term: Strong 



theoretical and phenomenological arguments indicate indeed that gluons do not bring any 
contribution of self-energy 66] . This term consequently vanishes in glueballs. For what con- 
cerns the instanton contribution, we can mention a previous study stating that instanton- 
induced interactions could significantly contribute in the ±+ glueballs [67|. But, to our 
knowledge, the inclusion of such interactions has not been elucidated yet within the frame- 
work of constituent gluon models and consequently, we will not consider any instanton-like 
contribution in glueballs. 

Finally, the basic remaining question to answer to is: Can glueballs (at least the low-lying 
ones) be described by a bound state of two constituent gluons? Roughly speaking, we gave 
an affirmative answer in Ref. Ssj]. We showed in this last work that the lattice QCD mass 
and wave function of the lowest lying ++ glueball (69] are compatible with those of a bound 
state of two constituent massless particles, interacting via a Cornell potential. The results 
of Ref. 68| are plotted in Fig. El The optimal potential, obtained by using the optimal 
values of the lattice QCD data, clearly exhibits a confining long-range part, and a rapidly 
decreasing short-range part. The errors on these lattice data allow the "true" potential to be 
located between two extremal curves, in the gray area. The curve we get in Ref. [68j can be 
compared to the potential for the ++ glueball in the present flux tube model. This state is 
a i = S = one. Consequently, only the contact interaction and the retardation term have 
to be added to the convoluted linear-plus- Coulomb potential. Using the approximation §§§ 
for massless particles, we find that the potential is of the form 



V ++(r) = a g r -3a s (r) 



U(r) + 



e -r/ Jg 



3a fl 

8/i 



(22) 
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In this potential, the tension a g for the flux tube in a glueball is related to the corresponding 
one in a meson, simply denoted as a, through a scaling law a g = Ca. The constant C is 
generally assumed to be given by 9/4 (Casimir scaling) [70j, or by 3/2 (square root of Casimir 



scaling) 



7l| . We assume here the Casimir scaling hypothesis, w\ 



70 



rich has been confirmed by 
22| ■ The choice C = 3/2 



lattice QCD calculations as well as by effective models [30, 
actually emerges from bag models-inspired approaches, that are not considered in this work. 
It is worth noting that a 3 factor is now present before ats(r), which corresponds to the value 
of the color operator in the one-gluon-exchange diagram between two gluons. Moreover, we 
showed in Ref. [73J] that, in the ground state, 



37ra g 



(23) 



The potential defined by relations (jzzj) and f[23]) can be compared to the results coming 
from lattice QCD: For «o = 0.400 and 7 9 = 7z = 0.940 GeV -1 , it is maximally located 
in the gray area, as it can be observed in Fig. El Interestingly, the effective gluon size is 
compatible with the one of massless quarks. Assuming the Casimir scaling hypothesis, the 
same parameters than for light mesons can thus be used. 



B. Mass spectrum 



The lowest-lying glueball states are given in Table IHIl and plotted in Fig. Our results 



75 



76l . [77J, excepted for the 2" 



+ 



are always located in the error bars of lattice QCD 
states. However, more states are present in our spectrum than in lattice QCD, i.e. J = 1 
states. Lattice calculations seem indeed to rule out the presence of 1 _+ and 1 ++ states 
below 4 GeV This can be qualitatively understood in terms of interpolating operators 
of minimal dimension, which can create glueball states, with the expectation that higher 
dimensional operators create higher mass states: The lowest states ++ , 2 ++ , h and 
2~ + are produced by dimension— 4 operators, while 1 ++ and 1 _+ states are respectively 
produced by dimension— 5 and dimension— 6 operators ItoT ]. Nevertheless, we can observe 
in Table llHl that our model predicts the existence of 1 _+ and 1 ++ states around 3 GeV. To 
our knowledge, no such experimental candidate has been found yet. The presence of these 
states in our model may actually be due to the use of spin degrees of freedom instead of 
helicity ones, which could be the correct ones if the constituent gluons are massless. Let 
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us note that this point is not alwa ys accepted in the literature: Many approaches involve 



indeed massive valence gluons 



27 



28 



29J. An intuitive argument is the following. If the 



spin is taken as internal degree of freedom, the Pauli principle states that S has to be even 
(0 or 2) for even £, and odd (only 1) for odd t. But, the combination of two gluons with 
the helicity h as degree of freedom gives the states h = A , s and 2, where O 5 ^ is a 
(anti) symmetric configuration. For even £, the possible helicities are s or 2, as for the spin. 
But, for odd £, h = A is the only possibility. The differences between spin and helicity 
degrees of freedom will thus mainly manifest for odd £. Interestingly, the 2 h states, which 
are missed in our model, are £ = 1 states. We thus think that, although the rough picture 
as a two-gluon bound state interacting via a Cornell potential seems to be valid, a more 
careful study of the relevant degrees of freedom for the constituent gluons is the key to 
a possible fully quantitative application of potential models to glueballs. The presence of 
spurious states around 3 GeV in our model could indicate what are the limits of the present 
approach. 

From the previous discussion, we can guess that the £ = glueball states are correctly 
described within our model, whether spin or helicity have to be used. These states are the 
++ and 2 ++ ones, which are compatible with the /q(1710) and the /2(2340) respectively. 
Consequently, our results support the identification of the / (1710) with a pure glueball 



state, in agreement with Ref. 
for the tensor glueball. 



611 ] . The / 2 (2340) should also be seen as a serious candidate 



VI. SCALAR GLUEBALL AND MESON MIXING 



In order to understand the numerous lig ht ++ states which have been experimentally 
detected, it was firstly suggested in Ref. 78J that these scalar states could be interpreted as 
mixings between pure \gg), \nn) and \ss) states. Many works have since been devoted to the 

sab, 



61 



79 



80 



81 



82 



study of the mixing between scalar mesons and glueball (see Refs. 
with different results following the approaches. In our framework, the masses of the pure 
++ states are 



M gg = 1.705 GeV, M sS = 1.528 GeV, M nn = 1.351 GeV. 



(24) 
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The structure of mixed states can be investigated thanks to a mixing matrix [79], which, 
following the lattice study of Ref . [82] , is of the form 



M 



99 

A M, 



A V2r A 


\y/2rX M m 
0.064 ± 0.013 GeV and r = 1.198 ± 0.072 



(25) 



In this last study, the mixing 



with A 

between pure states was used to understand the structure of the /o(1370), /o(1500), and 
/o(1710) mesons. However, following our previous conclusions, we assume that the /o(1710) 
is the pure glueball state, and that the /o(1370) is the pure nn state. 

We would like to add in the game the recently observed /o(1810) [84] . which can be con- 



sidered also as a good glueball candidate 



85] . It is worth mentioning that this interpretation 



is not unanimous: In some studies, the /o(1810) is interpreted as a mainly ss hybrid meson 
state 86] . However, hybrid mesons are not included in the present discussion. As we only 
consider the mixing between usual mesons and glueball, it is thus interesting to determine 
whether the / (1810), whose mass is 1.812 ± 0.044 GeV [84J], and the / (1500), whose mass 



is 1.507 ± 0.005 GeV 



48] , can be reproduced by the simple mixing matrix 



To do 



that, we begin by fixing the ratio r to the value 1.198, as it was com put ed by lattice QCD 

j, we will allow the 



calculations. Since we are interested in different states that in Ref. 
parameter A to be slightly different that the value 0.064 GeV. Interestingly, for A = 0.104 
GeV (35% higher than the maximal value of Ref [82j), we find three eigenstates defined as 
follows: 



Mi = 1.811 GeV, |V>i) = -0.883 \gg) - 0.325 \ss) - 0.338 \nn) 
M 2 = 1.502 GeV, \ip 2 ) = +0.237 \gg) - 0.931 \ss) + 0.277 \nn) 
M 3 = 1.272 GeV, |V> 3 ) = -0.405 \gg) + 0.164 \ss) + 0.899 \nn) 



(26) 
(27) 
(28) 



We see that, for this particular value of A, Mi is compatible with the mass of the /o(1810), 
which could thus be a mainly glueball state. Moreover, Mi can be identified with the 
/o(1500), with a dominant ss component. An additional light ++ meson is now present, 
which is still compatible with the /o(1370), because of its large error bar. Actually, the mass 
of this state have been measured in a rather large range, from 1.2 to 1.5 GeV [48]. Only 
keeping in mind the latest data, the pure nn state we previously found is merely compatible 
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with the state of Ref. [87|, which was detected at a mass of 1.350 ± 0.050 GeV. That clearly 



corresponds to the / (1370). But, in Ref. 



a state was detected with a mass equal to 



1.265 ± 0.065 GeV. This state was also identified with the /o(1370), but its mass is rather 
compatible with our M3 state. Our approach then suggests that two states hide behind 
the current /o(1370). The heaviest one is compatible with a mass of 1.370 GeV, while the 
lightest one has a mass around 1.270 GeV. Further experimental studies should confirm this 
hypothesis. 

It is worth mentioning that, for r = 1.270 and A = 0.077 GeV, which are the largest 
values allowed following the lattice study of Ref. [82|, one finds 

M x = 1.774 GeV, = -0.911 \gg) - 0.285 \ss) - 0.298 \nn) , (29) 
M 2 = 1.509 GeV, ^2) = +0.231 \gg) - 0.952 \ss) + 0.202 \nn) , (30) 
M 3 = 1.300 GeV, \ip 3 ) = -0.341 \gg) + 0.115 \ss) + 0.933 \nn) . (31) 

The previous conclusions remain qualitatively correct without fitting the parameter A, even 
if the agreement with the experimental masses is less accurate. 

VII. CONCLUSIONS 

We presented in this work a generalized version of the relativistic flux tube model. Not 
only the dynamical contribution of the flux tube were taken into account, but also the 
retardation effects and the spin-orbit interactions generated by the flux tube itself. We 
also added short-range interactions emerging from one-gluon-exchange process, and for the 
meson case, quark- self-energy and instanton-induced forces. All these physical ingredients 
allow to drop the inclusion of an arbitrary constant, which is usually used in potential models 
to fit the absolute mass scale of the mass spectrum. Moreover, all the parameters appearing 
in our model can be fitted on lattice QCD results, which strongly reduces its arbitrariness. 

We firstly computed the heavy meson spectra. The bb states are globally correctly repro- 
duced, as well as the cc states. In particular, our results suggest to identify the three new 
states V(3940), V(3940), and Z(3930), with the triplet x c j(2P). Let us note that a previous 



study already suggested that V(3940) and Z(3940) are the Xa{2P) and Xc2^P) |50j. 

Since the relativistic flux tube model is well defined for light quarks, and because the re- 
sults in the heavy meson sector were encouraging, we applied it to the light meson case. The 
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nn spectrum is rather well described, excepted that our model misses the isospin splittings, 
which we do not take into account in this approach. We compute that a ++ state is present 



at a mass of 1.351 GeV, in agreement with the /o(1370) [871 ] . Following our results, the 
/2(1950) also seems to be a pure nn state. As we did not include the possibility of mixing 
between nn and ss states, only three states, which can be seen as nearly pure ss, can be 
described: the 0(1020), the / 2 (2010), and the / 2 (2300). 

Finally, we tried to extend the model to glueballs. Assuming the Casimir scaling hypoth- 
esis, we got a spectrum which is mostly located in the error bars of the lattice QCD one. 
However, spurious states appear in our spectrum around 3 GeV, which are not present in 
lattice QCD. It consequently appears that the application of potential models to glueballs 
requires a particular care. Indeed, the gluons are massless, and consequently, their inter- 
nal degree of freedom should be helicity. However, the relativistic corrections we use are 
computed as if a gluon was a massive spin-1 particle, whose mass is the dynamical mass 
u = (p 2 )- We then obtain expressions involving spin degrees of freedom. As we argued 
in this work, a careful study of the states obtained either with helicity or spin should lead 
to a better understanding of the domain of validity of potential models. We leave such an 
analysis for future works. 

Our results concerning glueballs can fortunately be assumed to be valid for the lowest 
lying ++ and 2 ++ states. Our model then suggests that the / (1710) is a pure glueball and 
that the ^2(2340) is a very good candidate for the tensor glueball. Moreover, if we add the 
possibility of a mixing between scalar mesons and glueball, we find that our model is coherent 
with the /o(1810) and the /o(1500) as dominantly glueball and ss states respectively. We 
also predict a third state around 1.270 GeV, that could be a new scalar state, mostly nn, 
probably already detected [3] but wrongly identified with the / (1370). 
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a (GeV 2 ) 


0.185 


m b 


(GeV) 4.785 
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m c 


(GeV) 1.445 


lh (GeV" 1 ) 
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m s 


(GeV) 0.202 


7l = lg (GeV- 1 ) 


0.940 


m n 


0.000 


/ 


3.000 


m g 


0.000 


5 (GeV) 


1.000 







TABLE I: Numerical values of the parameters which are involved in our computations. m n com- 
monly denotes the mass of the n quark, which commonly denotes u or d quarks. 



> 
a 




FIG. 1: Effective quark-antiquark potential computed from lattice QCD data (full circles) [52|. 
Potential V h {r), given by Eq. ([19]), is plotted with a = 0.185 GeV 2 , a = 0.400, and -y h = 0.200 
GeV" 1 (solid line). 



29 



qq (n + l) 25+1 Lj M AM rft AM ret AM so AM qse (v og e{7) + ±a S {r)U(r)) 



bb 


l 3 ^ 


9.572 





-0.127 





-0.002 


0.018 






10.076 





-0.086 


o 


-0.002 


0.008 




l 3 Pn 


9.939 


-0.002 


-0.098 


0.005 


-0.002 


-0.028 




l 3 L>i 


10.180 


-0.003 


-0.084 


0.005 


-0.002 


-0.016 


cc 


l 3 Si 


3.240 





-0.166 


o 


-0.031 


0.053 




2 3 S'i 


3.814 


o 


-0.115 


o 


-0.029 


0.037 




l 3 Pn 


3.645 


-0.009 


-0.131 


0.027 


-0.030 


-0.080 




l 3 _Di 


3.937 


-0.017 


-0.112 


0.028 


-0.030 


-0.042 


nfi 


l 3 5i 


1.404 





-0.195 





-0.533 


0.075 




2 3 5'i 


2.067 





-0.130 





-0.352 


0.034 




1 3 P 
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TABLE II: Numerical evaluation of the different terms appearing in the mass formula (11 1H for some 
qq states. Parameters of Table U are used. The nonzero values are given in GeV. 
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FIG. 2: bb meson spectrum computed with our model and the parameters of Table[I] (empty circles) 



Our results are compared to the experimental data taken from the Particle Data Group 



48]. 
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FIG. 3: Same as Fig. [21 but for the cc mesons. More details about the X, Y, Z states can be found 
in Ref. 50] for example. Lattice results (crosses) are taken from Ref. [5 
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FIG. 4: Same as Fig. O but for the nn mesons. The instanton contribution (j20l) has been added 
to the 71 and states. 
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FIG. 5: Theoretical predictions for the masses of some highly excited nn mesons (empty circles). 
Parameters of Table [J were used. 
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FIG. 6: 
glueball 



ffective gluon-gluon potential computed from lattice QCD data [69] concerning the ++ 



681 ] . The optimal potential is plotted with a dotted line, but the error bars on the lattice 



results allow the true effective potential to be located in the gray area. Potential V ++(r), given by 
Eq. p2J), is plotted with a g = (9/4) 0.185 GeV 2 , a = 0.400, and lg = 0.940 GeV" 1 (solid line). 
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FIG. 7: Graphical representation of the results presented in Table IIIH The spectrum obtained 
with our model (empty circles) is compared to lattice QCD computations (crosses) and to possible 
experimental candidates mentioned by the Particle Data Group (full circles) [481 ] . 
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(n + 1) 2S + 




Model 


Lattice 


Experiment 


I 1 So 


0++ 


1.705 


1.710 ±0.050 ±0.080 [74] 


1.718 ±0.006 [48J 


1 5 S 2 


2++ 
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3.354 


3.040 ±0.040 ±0.150 [74] 




2 3 P 


o-+ 


3.449 


3.640 ±0.060 ±0.180 [75] 




1 X G 4 


4++ 


3.858 


3.650 ±0.060 ±0.180 [76] 
3.608 ±0.110 [77] 




1 5 G 3 


3++ 


3.895 


3.670 ±0.050 ±0.180 [74] 




1 3 F 2 


2-+ 


3.575 


3.890 ±0.040 ±0.190 [75] 




1^6 


6++ 


4.552 


4.339 ± 0.261 [77] 





TABLE III: Masses of some low-lying two-gluon glueballs computed with our model (third column) 
and the parameters of Table HI The quantum numbers of a particular J PC state are summed up in 
spectroscopic notation in the first column. Our results are compared to lattice QCD calculations 
(fourth column), and corresponding experimental candidates are suggested in the last column. All 
the masses are given in GeV. 
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